Phase boundaries of nanodots and nanoripples over a range of collision cascades 
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The nonlinear continuum model proposed by Cuerno and Barabasi is the most successful and 
widely acceptable theoretical description of oblique incidence ion sputtered surfaces to date and is 
quite robust in its predictions of the time evolution and scaling of interfaces driven by ion bom- 
bardment. However, this theory has thus far predicted only ripple topographies and rough surfaces 
for short and large scales, respectively. As a result, its application to the interpretation and study 
of nanodots, predicted by Monte Carlo simulations for, and observed in experiments of, oblique 
incidence sputtering is still unclear and, hence, an open problem. In this paper, we provide a new 
insight to the theory, within the same length scale, that explains nanodot formation on off-normal 
incidence sputtered surfaces, among others, and propose ways of observing the predicted topogra- 
phies of the MC simulations, as well as possible control of the size of the nanodots, in the framework 
of the continuum theory. 



For some time now, the scientific community has been 
captivated by possible cooperative behaviour exibited by 
surfaces sputtered by energetic ions, in which patterns 
are not only formed in the process of random ejection of 
surface particles, but that the orientation of the patterns 
can be tuned by varying the experimental sputtering pa- 
rameters 0-01 . Coupled with their nature as nanostruc- 
tures, these nanopatterns are of high technological im- 
portance due to obvious opto-electronic and other nano- 
device applications. Among the earliest theories put for- 
ward to explain the time evolution of ion eroded surfaces 
was the linear continuum theory of Bradley and Harper 
PJ, which describes the surface topography at any time 
instant as the result of an interplay or competition be- 
tween the destabilizing erosion process, that creates an 
instability in which troughs are eroded in preference to 
crests, and the stabilizing processes of surface diffusion. 
Although, this explains the short length scale properties 
of sputtered surfaces observed in the period, it does not 
explain the larger length scale roughening behaviour. 

According to Cuerno and Barabasi, the latter case is 
the result of nonlinear effects and the sputter noise in 
the sputtering process While this defines new scal- 
ing regimes absent from the linear theory, as well as rip- 
ple topographic regimes agreeing with the linear theory, 
a number of unresolved problems persist. For instance, 
the scaling properties of the different scaling regimes are 
unknown. Furthermore, the nonlinear theory, like the 
linear theory, has thus far predicted mainly ripple to- 
pography at short length scales. Also, calculations are 
usually performed for isotropic distribution of the energy 
of the impinging ion, mainly for ease of exposition of the 
theory, and a few cases of anisotropic energy distribution. 
These are inconclusive for certain important cases as we 
shall see below. 

Meanwhile, a number of experiments of normal inci- 
dence sputtering have demonstrated the existence of a 
different kind of surface morphology, dot topographies, 



on semiconductor and amorphous surfaces at nanometre 
length scales with the presence of a characteristic length 
scale in the system [fj-|2|. Moreover, recent Monte Carlo 
simulations [13, [H| reported the existence of these dot to- 
pographies for off-normal incidence sputtering of amor- 
phous substrates at collision cascade parameters different 
from the existing continuum theory calculations. 

In Ref. [Tot six different topographic regions were re- 
ported as to be expected for early times in the sputter- 
ing process. When considering the topographies at later 
times as well as the nature of the nano structures, this 
six regions reduce to three, through a merger of the first 
four (I, II, III, IV). 

Although, a number of targeted theoretical descrip- 
tions focussed at accounting for nanodot characteristics 
have been proposed, 12Hl4j the actual formation of these 



nanodots from oblique incidence sputtering remain un- 
clear and is not yet understood, since the continuum the- 
ory has so far not predicted anything other than ripple to- 
pographies for off-normal incidence. This is still an open 
problem, which we investigate in this paper by provid- 
ing phase diagram calculations of the continuum theory 
yet unreported for anisotropic distribution of the energy 
of the impinging ion. In particular, we provide those 
necessary to resolve the unanswered questions about the 
prediction of the continuum theory as regards non-ripple 
morphologies. 

We propose a different interpretation of the continuum 
theory which, among other explanations, indicates that 
for anisotropic distribution of the energy of the imping- 
ing ions, the presence of a characteristic length scale may 
not predict the separation of ripple crests or troughs but 
that of dots. We discuss the accessibility of any of the 
regions studied to experimental probing; propose an ex- 
planation for the formation of the dots that arise from 
oblique incidence ion sputtering on the basis of the con- 
tinuum theory, and a possible way of achieving or con- 
trolling required dot size. 
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The continuum theoretical description of interface 
morphology in terms of deterministic and stochastic par- 
tial differential equations is a powerful and successful tool 
for understanding the behaviour of diverse interface phe- 
nomena. For the specific case of ion sputtered surfaces, 
the distribution E(x) of the energy E of the incident ion 
to a surface particle located at position x = (xi, X2, X3) is 
assumed in the continuum theory to be of the Gaussian 
form fil: 
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where a and p are the widths of the distribution paral- 
lel and perpendicular to the ion beam direction, respec- 
tively. The erosion velocity v oc dth, by definition, fol- 
lowing which the dynamic evolution of the surface height, 
h(r,t), at nanometre length scales is, for most cases, gov- 
erned by a Kuramoto-Sivashinsky type stochastic partial 
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differential equation [5| 

d t h(r, t) = -v + (d x h(r, t) + s x d xx h(r, t) + q y dy y h{r, t) 
+ Vx [d x h(r, t)] 2 + r} y [d y h{r, t)f - DV 4 h{r, t) + f3. (2) 

Vq is the erosion velocity of a flat surface, £ is a propor- 
tionality constant related to the local surface slope along 
the x-direction, v x and v y are the (linear) surface tension 
coefficients, A^ and Aj, are the nonlinear coefficients, D is 
the surface diffusion coefficient, and /3 is the normal noise 
term, which is assumed to have a Gaussian distribution 
with zero mean. 

Using the convenient notation, 
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where F is the ion flux and J is the proportionality con- 
stant between the power deposition and the rate of ero- 
sion, we provide the coefficients, for ease of reference, as 
follows. dEl: 
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Non-linear effects are believed to be irrelevant at nano 
scales, hence, at such length scales Eq. [5] predicts the 
presence of a characteristic length scale T = y/D/\<;\ in 
the system 0, [B| which manifests as periodic structures 
(e.g. in the separation of ripple crests/troughs); where 
\<;\ is the largest absolute value of the negative surface 
tension coefficients. Thus, if neither of \^ x \ and \<; y \ is less 
than zero, the characteristic length scale is absent. And, 
since the ripple wavelength is A = 2tt\/2T, no ripples are 
formed. In other words, the present continuum theory 
interpretation is that we either have ripples or not for 
oblique incidence, whereas reports of the possibility of 
other topographies have emerged. 

In order to clarify this we obtain the phase diagrams of 
the continuum theory for the asymmetric energy distri- 
bution cases a — — 5 and p = — 5, reported in Refs. [13 
The phase diagrams have been obtained from the 
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variations of S x ,Sy,r) x , % as functions of a and p (Figs. Q] 
and [5]); and as functions a and 9 (Figs. [31 andSJ). A cal- 
culation of the coefficients (as functions of a and p) have 



r 



revealed the same number (three) of topographic regions 
(Fig.[T]), as in the simulations at later times. Since these 
coefficients are also functions of a and 9 it is necessary to 
obtain the phase diagram in these regard as well. Indeed, 
as we show below, the phase diagram in terms of a and 
9 allow for a wide range of possibilities of the coefficients 
and, hence, a large number of topographic regions. The 
regions to be encountered below in any of the two cases 
are as defined in Table HI As will be seen in what follows, 
if we ignore the nonlinear coefficients then there are only 
three possible regions, two of which describe ripples with 
either of the two possible orientations and the remain- 
ing one describing the situation in which neither of the 
surface tension coefficients is negative. 

In the phase diagram presented in Fig. [1] obtained 
from the continuum theory, ripples are oriented along 
the y-axis in the three regions 2, 4, and 9. Whereas, in 
the simulation ripples are oriented perpendicular to the 
ion beam direction for the same sputtering parameters. 
Assuming that the projection of the ion beam direction 
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TABLE I. Definition of the regions found in the calculations, 
as described in the text. 



Region Relative Signs of the Coefficients 
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FIG. 1. Phase diagram for collision cascade parameters a and 
p ranging from to 5. 8 = 50°, a = 6.0. The three regions 2, 
4, and 9, are as defined in the text. 



in the simulation corresponds, in the reference frame of 
the continuum theory, to a straight line segment par- 
allel to the x-axis, we have an agreement between the 
continuum and the discrete theories that enables us to 
interprete the results of the continuum theory as regards 
dot topographies. For instance, since the disagreement 
here is in the topography of region 9, then the nonlinear 
coefficients must be dominant and we explain their role 
below. Note that the colouring in Fig. [T] (online version) 
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FIG. 2. Boundary shifts in the phase diagram, arising from 
varying 9, for the same parameters as in Fig.fl] Left: 6 — 70°; 
Right: 6 = 30°. 

is to further enhance the distinction between boundaries, 
and not coded for any specific region. 

Until now, phase diagram calculations have not been 
done for this oblique incidence region and the only contin- 
uum theory interpretation of the results of applications of 
Eq. [2] to an understanding of the time evolution of sput- 
tered surfaces is one that merges two relatively different 
length scales, one short (of the order of 1/im, the other 
long (of the order of tens of fim). This interpretation 
is, of course, very crucial to an understanding of transi- 
tions that do occur, has been experimentally verified, and 
are exhibited in the surface topography and the scaling 
behaviour reported by diverse experiments at different 
length scales. Here, we provide the results of these yet 
unreported calculations and on their basis propose an- 
other interpretation of the result of an application of Eq. 
[2]within the same length scale, which, we argue, accounts 
for the unexplained phenomenon of oblique incidence dot 
formation. 

Note that the values of the collision cascade parame- 
ters at the phase boundaries are a bit different to those 
in the simulation results. To investigate this we perform 
calculations for the phase diagram at different 9 [see Fig. 
[2], and observe shifts in these boundaries, which indicate 
the possiblity of a quantitative agreement with the val- 
ues reported in the simulation. The general result for the 
three regions of Fig.[lJ when considering both a and 0, is 
shown in Figs. [3] and 0] (a) and (b), for a — 3.3, p = 4.5; 
a = 1.6, p = 3.3; and a = 2.0, p = 1.0, respectively. A 
transition from region 2 to 4 of Fig. [U when considering 
a fixed a and varying p, is due to the change of sign of 
the nonlinear coefficient r/ x , associated to increasing local 
surface slopes along the x-axis, which is a precursor to 
the change from region 4 to 9, in which case there is no 
directional change in nonlinearity but only a change of 
sign of the surface tension coefficient q x . This is impor- 
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FIG. 3. Phase diagram for the isotropic case a = 3.3 and 
p — 4.5, representative of region 9 of Fig. [1] for varying 6 
and a; ^(horizontal axis) ranging from to 90, and a (vertical 
axis) ranging from to 6. The regions found are as labelled 
in the figure, and defined in the text. 



tant as it provides an explanation of dot formation. 

A negative surface tension coefficient is representative 
of the instability arising from the sputtering process in 
which troughs are eroded in preference to crests. A posi- 
tive surface tension coefficient would then imply a neglect 
of troughs in the erosion process. On the other hand, 
a negative nonlinear coefficient implies that the height 
evolution increases as local surface slopes increase, and 
vice-versa. This means that the interplay that leads to 
ripple formation is enhanced or countered depending on 
the relative signs of the nonlinear coefficients which is 
capable of creating a further instability that disturbs the 
interplay. Following this, pattern formation in the 12 re- 
gions highlighted above are as tabulated below in Table 
HI1 Note that due to the interplay between u x and v y 
periodic structures with either of two possible orienta- 
tions are always present, except if one of v x or v y is zero. 
Details of how to calculate the quantities (e.g. ripple am- 
plitude, growth rate, etc.) in Table UU can be found in 
Refs. [Etui. 

In particular, for an explanation of region 9, trough 
erosion along the x-axis is not favoured and the height 
evolution along the x-axis decreases. Since, the erosion is 
a stochastic process, and erosion along the x-axis, for this 



FIG. 4. Phase diagram for representative colision cascade 
parameters of the two remaining regions of Fig. [TJ and for 
varying 6 and a; ^(horizontal axis) ranging from to 90, and 
a (vertical axis) ranging from to 6. (a) a — 1.6 and p — 3.3; 
(b) a = 2.0 and p = 1.0. The regions found are as labelled in 
the figure, and defined in the text. 



region, is much reduced in comparison to that along the 
y-axis, the continuity of eroded troughs along the y-axis 
is broken and instead of long grooves we have pits inter- 
spersed with isolated protrusions which together make 
dot topography. 

Based on this explanation, then dot size depends on the 
relative magnitudes of the nonlinear coefficients, which 
again is dependent on the sputtering parameters. Thus, 
for a preferred dot size and growth with time, one would 
need to strike the right balance between the appropriate 
choices of material (which influences a, a, p, etc) and 
sputtering conditions such as ion incidence, temperature, 
etc, according to the phase diagram of Figs. [TJ [2j [3l and|4] 
There are a few cases with ? x > 0, q y < 0, rj x < 0, rjy = 
in Figs. [5]and[31 <; x > 0, <; y > 0, 7] x = r\ y = in Fig. [3] and 
a few case of region 12 in Figure ??. These do not appear 
in the phase diagram because they are exceptional cases 
which are too tiny to be noticed. 

On the accessibility of some of the regions studied here 
to probing by experiments, we recall the assumption of 
the continuum and discrete 1(1 11, 12, HI theories 
that the impinging ion penetrates a distance a into the 
material, after which it distributes its energy according 
to a Gaussian distribution, which is a simple represen- 
tation of the geometry defined by the collision cascades 
triggered by the impinging ion. Thus, it is possible to 
design the shape of this geometry through controlled de- 
fect creation or doping in the first few surface layers of 
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TABLE II. Definition of the regions found in the calculations, as described in the text. 
Region Topography 

1 ripples oriented along x with much less prominent underlying periodic structure along y 

2 ripples oriented along y with much less prominent underlying periodic structure along x 

3 ripples oriented along x with shorter amplitude and (possibly prominent, depending on the relative size of ? y ) 
periodic structure along y 

4 almost pure ripples (i.e. very little or no sign of an underlying periodic structure) oriented along y with normal 
amplitude growth 

5 almost pure ripples oriented along x with normal amplitude growth 

6 ripples oriented along y with shorter amplitude and (possibly prominent, depending on the relative size of ? T ) 
periodic structure along x 

7 possibly a rough surface, or low amplitude ripples oriented along x; depending on the relative strengths of the 
competing factors 

8 ripples oriented along y with slow or no amplitude growth; or rough surface if the nonlinearities cancel out or soften 
the surface tension instability 

9 dots with underlying periodic structure oriented along y 

10 short ripples oriented along y, or dots of lower growth rate and underlying structure oriented along y; depending or 
the relative strengths of the competing factors 

11 ripples oriented along x with slow or no amplitude growth; or rough surface if the nonlinearities cancel out or 
soften the surface tension instability 

12 dots with less prominent or no underlying periodic structure 



width roughly about the penetration depth a. 

In summary, we have proposed a new interpretation of 
the continuum theory that explains surface topographies 
yet unaccounted for and also agrees with the previous 
conceptual framework of the theory. This new insight 
considers the stochastic time evolution equation within 
the same lengthscale that is of general interest in exper- 
imental studies and characterization of the nanostruc- 
tures. A possible means of achieving the sputtering pa- 
rameters neccessary for experimental observation of the 
topographic regions of present interest, e.g. the nanodot 
regions, is discussed. 

The author thanks Alexander Hartmann and Reiner 
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